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I The slope of the FP therefore departs from the virial expectation of oc (To(S)^ 

' at all optical and near-infrared wavelengths, which could be a result of the variation 

• of M/L along the elliptical galaxy sequence, or a systematic breakdown of homology 
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1. Introduction 



Correlations among the properties of elliptical galaxies have been used both as measures of 
their homogeneity as a population and as indicators of the distances of individual galaxies. The 
discovery of a color-magnitude effect ( |Baum 1959 ) was followed by the measurement of relative 
distances of galaxies and clusters using the color-magnitude relation ( ^andage 1972 ; Visvanathan 



& Sandage 1977; Sandage & Visvanathan 1978a, b). The relation's small scatter was used as a 
constraint on elliptical galaxy formation time-scales ([Bower, Lucey, &: Ellis 1992| ). The correlation 
between luminosity and velocity dispersion ( Faber Sz Jackson 19761) was likewise used as a distance 
indicator ([Tonry fc Davis 1981 ; Dressier 1984| ), but Terlevich et al. (1981) discovered a weak 
correlation between the residuals of the relation and Mg2, implying that there might be a second 
parameter which contributes to the intrinsic scatter of the Faber-Jackson relation.^ Subsequent 
studies (Dressier et al. 1987; Djorgovski & Davis 1987) that included large samples of elliptical 
galaxies found that surface brightness was a second parameter which caused a large portion of 
the scatter in the Faber-Jackson relation. In this perspective, the intrinsic properties of elliptical 
galaxies are only found to lie on a plane within the three-dimensional parameter space of the 
observables. This Fundamental Plane (FP) is thus a set of bivariate correlations between the 
observed properties of elliptical galaxies; the color-magnitude and Faber-Jackson relations are 
projections of that plane onto two of the three axes. 

The importance of the exact form of the slope of the FP was immediately identified as possibly 
providing a strong constraint on the mass-to-light ratios (M/L) of elliptical galaxies ( Dressier 
al. 1987; pjorgovski fc Davis 1987 ). In particular, the virial theorem produces a prediction that 
Teff oc cro(S)gg^ (where r^s is the effective radius, ctq is the central velocity dispersion, and (S)eff 
is the mean surface brightness enclosed within res) if two assumptions are made: (1) M/L is the 
same for all elliptical galaxies, and (2) elliptical galaxies form a homologous family in their scaling 
properties. The latter assumption was generally taken to be true, hence virtually all researchers in 
the last decade have proceeded to explore the effects of the variations of M/L implied by the FP 
correlations. For example, the slope of Dn oc ctq'^ from Lynden-Bell et al. (1988) in the -B-band 
implies that M/L varies systematically with the galaxy's luminosity as M/Lb oc L^^"^. 

More recently, there have been a variety of theoretical (Capelato, de Carvalho, & Carlberg 
1995, 1997; Ciotti, Lanzoni, &: Renzini 19961) and observational ( praham fc Colless 1997 ; Busarello 



et al. 1997) investigations into the question of whether or not elliptical galaxies form a homologous 
scaling family. The results of these studies are not yet clear, but they seem to imply that structural 
deviations of the light profiles of ellipticals from a homologous family cannot affect the FP 



''Dressier et al. (1987) describe how this evidence for the bivariate nature of elhptical galaxies found by Terlevich 
et al. (1981) might actually have been driven more by distance errors in the Terlevich et al. sample, which was drawn 
primarily from the field, and from a surface brightness correlation which was later found to be the second parameter. 
After correcting for surface brightness effects, Dressier et al. found little or no correlation among the residuals of the 
Faber-Jackson relation and Mg2. 
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appreciably ( [Graham fc Colless 199?] ), while dynamical deviations from a homologous family can 
(Capelato et al. 1995; Busarello et al. 1997; cf. Graham fc Colless 1997| ). Underlying all of these 



studies is an important point: if there are significant and systematic deviations from a homology, 
then these deviations should be strictly independent of wavelength observed in constructing the 
global photometric parameters which enter the FP. 

In a broad sense, some of these possible implications of the FP make specific predictions which 
can be tested by obtaining additional data. For example, if the form of the FP is a direct result of 
a dependence of M/L on L due to variations in the stellar populations parameters of age and/or 
metallicity, then observations in the near-infrared should show a significantly different form for 
the FP correlations as the 2.2//m light is far less sensitive than optical light to line-blanketing and 
somewhat less sensitive to age effects. Alternatively, if the origin of the FP is due to a systematic 
deviation of elliptical galaxies from a homologous family, then the exact form of the FP should 
be independent of wavelength. It is also possible that some combination of these effects could be 
required by a simultaneous analysis of the FP at optical and near-infrared wavelengths. 

The present paper is an attempt to address these possible origins for the FP correlations by 
exploring their form using near-infrared imaging data. Early work on the FP in the X-band 
was done by Recillas-Cruz et al. (1990, 1991), who obtained aperture photometry for galaxies in 
Virgo and Coma, and Djorgovski & Santiago (1993), who used aperture photometry from Persson, 
Frogel, & Aaronson (1979). Both studies relied on optical estimates of r^^. It follows directly from 
the existence of color gradients that Teg should be smaller for longer wavelengths and this point 
will be shown explicitly in a future paper ( pahre, de Carvalho, fc Djorgovski 1998 ) for these optical 



and near-infrared data. The present paper is more than just a revisiting of the X-band FP — it 
is an attempt to study the global properties of elliptical galaxies using near-infrared photometry 
that is fully independent of the optical photometry, while at the same time using a methodology 
for deriving global photometric parameters that is identical to the method at optical wavelengths. 

An imaging survey of this kind and scale has only recently become possible with the 
introduction of large format IR detectors (256 x 256 pixel^). This project was initiated during the 
commissioning phase of a wide-field, near-infrared imaging camera ( [Murphy et al. 1995| ) for the 
Palomar 60-inch Telescope. In the first paper of this series, early results from this survey on the 
i^-band FP ( Pahre, Djorgovski, fc de Carvalho 1995[ ) indicated that there is a modest change in 



slope from the optical to the near-infrared, but not nearly as much variation as would be expected 



if stellar-populations alone were the cause of the slope of the FP ( Pahre fc Djorgovski 1997 



Pahre, Djorgovski, fc de Carvalho 1997 ). The full K-hand survey and the complete catalogs of 
global properties are described in the second paper of this series ( pahre 1998[ ) and are summarized 
in §^. All of the data contained in the previous contributions (Pahre et al. 1995, 1997; Pahre & 
Djorgovski 1997) were re-calibrated, some were re-reduced, and the global photometric parameters 
re-derived in a consistent manner as described in Pahre (1998). The FP correlations and their 
many projections are derived in §H as a way of exploring various aspects of these near-infrared 



data. The Mg2-a"o relation is constructed for the same galaxies in §3.4. Simple models will be 
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compared to these results in but it will be shown that naive models cannot explain the many 
observed properties of ellipticals nor are such models unique. 



2. Description of the Data 

The data used for this paper derive from Pahre (1998). That paper presented near-infrared 
K-band imaging of 341 early-type galaxies, and used those data to measure the global photometric 
parameters of the half-light effective radius reg, the mean surface brightness (//)efr enclosed by that 
radius, the total magnitude Ktot, and the diameter Dk at which the mean surface brightness, fully 
corrected for cosmological effects and extinction, drops to 16.6 mag arcsec"^. The latter quantity 
is an analog of the 5-band -D„ parameter defined by Dressier et al. (1987). The near-infrared 
data were corrected for the effects of seeing. As shown in that paper, the random uncertainties of 
the measured quantities are: 0.06 dex on Teg; 0.21 mag on {fJ-)eS', 0.09 mag on i^tot! 0.010 dex on 
Dk', and 0.015 dex on r^s — 0.32(/x)efr, the quantity which will enter the FP. 

The galaxies in that sample are primarily drawn from nearby rich clusters of galaxies, 
although additional galaxies were added from groups and the general field. The galaxies were not 
selected according to any explicit criteria of completeness (such as a magnitude-limited sample 
would be), but by the availability of companion optical imaging and spectroscopic data. The 
primary goal of this effort was to provide a large sample of galaxies for which the variations of the 
FP correlations between the optical and near-infrared wavelengths could be explored. The data 
probe the full range of properties [ves, Ktot, ctq) displayed by the family of giant elliptical galaxies, 
and a significant portion of the sample is comprised of SO galaxies. 

The optical global photometric parameters {r^g, {^J')e{i, and Dn) and spectroscopic parameters 
(central velocity dispersion do and Magnesium line index Mg2) were compiled from the literature 
by Pahre (1998). All of the photometry were drawn from the U, B, g, V, Rq (or r), or Iq 
bandpasses and converted to V magnitudes for a general catalog. Furthermore, separate catalogs 
were constructed for individual comparisons to preserve the wavelength information for each 
optical bandpass. The spectroscopic parameters were corrected for observed aperture size effects 
to a common physical scale of 1.53h^^ kpc; small offsets between data sets have been applied 
according to prescriptions developed by other authors in the literature. The values were then 
averaged to reduce the random uncertainties and minimize systematic errors due to some data 
sets. Of the 341 galaxies imaged in the K-hand, 95% have velocity dispersions, 69% have Mg2 
indices, and 91% have optical photometric parameters (either r^s or -D„). The typical uncertainties 
for individual measurements of do and Mg2 are 0.04 dex and 0.013 mag ( Smith et al. 1997|) , 



respectively. A substantial fraction of the entire sample has more than one measurement of these 
parameters which were then averaged, so these two uncertainty estimates can be taken as a 
universal upper limit to the measurement uncertainties. 

Many, if not most, of the literature sources suggest that their velocity dispersions are less 
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reliable below 100 km s~^, but a bias in the slope of the FP can be introduced by imposing a 
cut on do. It will be important to investigate what effect changing the lower cutoff for (Tq has on 
the slope and scatter of the FP. Small measures of the effective radius, such as r^s < 1 arcsec 
(the median seeing was 1.34 arcsec FWHM), have large random uncertainties and probably 
substantial systematic errors, and should probably also be discarded. In the sample, 4% of the 
galaxies have morphological type SBO or later, and another 4% have SO/a type; caution should 
be exercised when studying the global properties of these galaxies. One galaxy (D45 in cluster 
Abell 194) appears to be a misidentification either in the optical or near-infrared as evidenced 
by its color {V — K) = 1.28 mag, which appears to be much too blue compared to the mean 
(y — K) = 3.15 mag for the entire sample. Five galaxies in the Virgo cluster were removed from 
the sample, as their accurate distances as derived by the surface brightness fluctuations method 
(provided by J. Blakeslee and J. Tonry; see Tonry et al. 1997| ) show that they are either in the 
background W Cloud (NGC 4168, NGC 4261, and NGC 4365) or the foreground (NGC 4660 and 
NGC 4697). 



3. Analysis of the Elliptical Galaxy Correlations 

3.1. The Near Infrared Fundamental Plane 

Galaxies were drawn from the sample described in §§. Only those galaxies residing in a 
cluster or group with four or more observed galaxies were included in the FP fits, resulting in 
16 clusters/groups and 249 galaxies with ctq > 1.8. Two of the five Leo I Group galaxies were 
excluded as a result of these selection criteria, although the remaining three galaxies were retained 
in the sample for completeness. 

The "standard" FP equation is usually written as 

log Teff (arcsec) = alogcro(km s^"*^) + 6(^)cfr(mag arcsec"^) + q (1) 

where a is usually identified as the "slope" of the FP and Cj are the "intercepts." The intercept 
of the relation will vary with distance. The galaxies in each of the 16 groups and clusters are 
assumed to lie at the same distance, hence there are i = 16 different intercepts c,. 

The Equation || was fit by minimizing the sum of the absolute deviations of the points 
orthogonally from the relation using the program GAUSSFIT ( |Jefferys et al. 1987 ). During the 



first iteration, two points which are outliers (D9 in Cen30 and S201 in Hydra) were identified and 
excluded from the analysis. The resulting FP in the near-infrared ET-band is 

log Teff (arcsec) = 1.531ogcro+ 0.314:{fiK)eS + Ci N = 251 rms = 0.096 dex 

±0.08 ±0.011 ^ ^ 



The uncertainties in the coefficients were determined by 100 iterations of bootstrap resampling of 
the data-points. The individual intercepts for the fit, and the rms about the fit for each cluster 
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or group, are listed in Table |l[ Since the rms is quoted in units of logrgfr, the uncertainty on each 
intercept is therefore rms/\/N — 1. The relation in Equation ^ is equivalent to the scaling relation 

^ ^ ^1.53±0.08/v \-0.79±0.03 

Changing the lower cutoff for do from 1.8 dex to 2.0 dex changes the value of a by < 0.01 dex, 
changes 6 by < 0.001 dex, reduces the scatter by 10% to 0.089 dex, and excludes 23 galaxies 
(9% of the total). Hence, the solution to the FP is robust to the lower ctq cutoff, although the 
galaxies with the lowest ctq appear to contribute the largest to the observed scatter. Minimizing 
the orthogonal variance (instead of the absolute value of the deviation) from the FP relation with 
(To > 1.8 dex results in a small change in the slope of the FP to: 

logreff(arcsec) = 1.631ogcro+ 0.320{fiK)eS + Ci N = 251 rms = 0.099 dex 

±0.06 ±0.008 ^ ^ 

Since the coefficients in Equations Q and ^ are equivalent within the uncertainties, the fit to the 
FP is insensitive to the exact fitting method. The method of minimizing the absolute value of the 
orthogonal deviation from the fit is to be preferred, however, as it is less sensitive to outliers. 

The simultaneous fit to all clusters is displayed in Figure |l| with the data subdivided into 
the 16 individual clusters or groups. It is clear from this figure that the simultaneous fit is a 
representative description of the properties of the early-type galaxies in all of the clusters. There 
is no clear deviation from this mean relation. 

The 11 clusters with more than ten galaxies were fit individually to Equation || as a test of 
the universality of the FP relation. The difficulty with all 11 separate fits, however, is that the 
number of galaxies in each cluster is small enough that the slope of the relation is not accurately 
determined. Once again, the uncertainties on the fitted coefficients a and b have been determined 
using the bootstrap procedure. In eight of the 11 cases, the fits have a within one standard 
deviation of the value a = 1.53 from the simultaneous fit, suggesting both that the uncertainties 
are reasonably estimated and that better fits are limited by the number of galaxies per cluster. 
These individual fits are listed in Table |l[ 

The same 11 clusters were also fit individually by constraining b = 0.314 from the previous 
simultaneous fit. This is possible because virtually every study of the FP (optical and near- 
infrared) obtains a similar value for this parameter, hence it should be possible to constrain its 
value a priori. These fits show significantly smaller uncertainty in their determination of a than 
the unconstrained fits, and are listed in Table |l|. In this case, seven of the 11 clusters have a slope 
a within one standard deviation of the value from the simultaneous fit. 

The adopted form of the FP in Equation |2| is plotted for all 301 galaxies in these 16 clusters 
and groups in Figure |2|, both in face-on and edge-on perspectives. While the edge-on view with 
log Res as the ordinate^ is the most common method of displaying the FP, the edge-on view 



^In this paper, a distinction will be made between the angular effective radius Tcb, measured in arcsec, and the 
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with log i?cfr — 0.314(/u)cfT as the ordinate is easier to interpret. Virtually every study of the FP 
(optical and near-infrared) obtains the same relationship between R^g and (/x)efr, but there may 
be significant variation in the relationship between (log-RgfT — 0.314(//)eff ) and logao, depending on 
wavelength. Furthermore, this edge-on perspective of the FP, seen from its short side, separates 
the correlated measurement errors in reff and (//)efr from the independent measurement errors in ao- 
The FP in physical units as plotted in Figure ^ is 

logiicfr(/i75^ kpc) = 1.53 log cro + 0.314(/x)cfr-8.30. 

In the face-on view of the FP in Figure ||(c), it is seen that galaxies do not uniformly 
populate this planar surface. While the -fC-band data in this paper are not drawn from a strictly 
magnitude-limited sample, they do behave as though a Ktot ^ 13 mag limit were imposed. Most 
of the galaxies are found to have 15 < {ijlk)cE < 18 mag arcsec"^ (long-dashed lines), although 
there are no clear selection effects causing this distribution of galaxy properties. Furthermore, 
there are no galaxies with properties in the upper-right portion of the figure, which could be 
caused by the lack of galaxies with central velocity dispersions (Jq > 400 km s^^, although there is 
no selection limit imposed on this portion of the FP. 



3.2. The DK-ao Relation 

Dressier et al. (1987) introduced a parameter D„ which was defined as the diameter at 
which the circular mean surface brightness (fully corrected for cosmological effects and Galactic 
extinction) dropped to a fiducial value. This parameter was chosen, in effect, to be a combination 
of the Tefj and (/u)eff terms in the FP correlations, thereby simplifying the FP to a Dn-ao relation. 
They defined this fiducial surface brightness to be 20.75 mag arcsec"^ in the i?-band: the surface 
brightness was faint enough that Db was typically much larger than the seeing disk, while bright 
enough that interpolation (rather than extrapolation) was used to evaluate Db from their aperture 
magnitude data.^ 

Lucey &: Carter (1988) defined an equivalent Dy parameter in the T^-band to be the diameter 
at which the mean surface brightness drops to 19.8 mag arcsec"^ (this assumes a mean galaxy 
color of {B — V)o = 0.95 mag). Smith et al. (1997) defined Dr for the i?c-band to be at 
(M)eff = 19.23 mag arcsec^^ (assuming {V — Rq) = 0.57 mag), J0rgensen et al. (1995a) defined 
the Gunn r-band Dr = 2r„ to be at (^)eff = 19.6 mag arcsec"^ (assuming (V — r) = 0.2 mag), 
and Pahre (1998) defined the -fC-band Dk to be at (/x)efr = 16.6 mag arcsec"^ (assuming 
(V — K) = 3.2 mag). By using typical colors for early-type galaxies in constructing these 
definitions, the average value of D measured for a sample of galaxies should be approximately 



effective radius R^s in physical scale, measured in kpc. 

*The notation ol Db will be adopted for the rest of the paper to distinguish the D„ parameter as defined in the 
_B-band from the equivalent diameter as defined in another bandpass. The name "D-cro relation" will refer to the 
correlation in all bandpasses. 
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independent of bandpass. The slope of the D-gq relation may vary between bandpasses, however, 
causing there to be a systematic variation of the D parameter as a function of a in different 
bandpasses (while keeping the mean D similar). For example, if Dy oc o""^ and oc cr'^^, then 
log Dy — log Dk = {ay — ax) log a + constant. 

The Dk-cq relation was fit for the galaxies in the same 16 clusters and groups as in § p.lj , 
excluding galaxies using similar criteria {Dk < 2 arcsec, logco < 1.8, late types), and using the 
bootstrap method to estimate uncertainties in the fitted coefficients. The best fitting relation is 

logDi^(/i75^ kpc) = 1.62 log CTo - 2.984 iV = 252 rms = 0.112 dex 

±0.07 ^ ^ 

and is displayed in Figure |^. The slope of this relation is consistent, given the estimated 
uncertainties, with the full FP relation in Equations ^ and |^. The scatter of the Dk-tq relation, 
however, is 15% higher than the i^-band FP, despite the fact that the measurement uncertainty 
of Dk is actually smaller than that for r^s — 0.32(//)efT ( Pahre 199^ ). This is most likely due to 



the fact that the Dk-ctq relation is nearly, but not quite, viewing the FP edge-on. 



3.3. The FP As Seen in k Space 

Since elliptical galaxies only populate a plane in the three-dimensional space of the 
observables {ReS, (/^)eff5 cro)) it is straightforward to construct an orthogonal transformation from 
this observed coordinate system to another one which might facilitate a physical interpretation 
of the underlying parameters. Since there are many possible transformations to accomplish this, 
the chief practical difficulty is identifying a transformation which permits a robust and unbiased 
physical interpretation of the true galaxy properties ( pjorgovski, de Carvalho, fc Han 198^ ). A 
transformation of this kind was proposed by Bender, Burstein, & Faber (1992): 

Ki = (21ogao + logi?efr)/V2 

K2 = (21ogao + 0.8(/i)efr-logi?cfr)/V6 (5) 
K3 = (21ogao + 0.4(/i)efr -logi?efr)/V3 

where the quantities Ki were constructed with the hope that ki would be proportional to mass, 
At3 would be proportional to mass-to-light ratio, and K2 (which is required to be orthogonal to 
Ki and K3) would be proportional to the product of mass-to-light ratio and the third power of 
mean surface brightness. This "K-space" is displayed in Figure |^ for the JC-band survey. Given 
the above interpretation of ki and K3, the fitted line between these two variables 



K3 = 0.147^1 + 5.721{fiK)eS N = 251 rms = 0.068 dex 
±0.011 ±0.038 



(6) 
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then implies that the "observed mass-to-hght ratio" ^ in the K-hand varies as 
(M/Lk) oc M°-i^'^=t°-°i^ oc L°^^72±o.oi3 ^ rpj^-g cQj^ciygiQj^ ig dependent on ehiptical galaxies 
forming a dynamically homologous family in which the central velocity dispersion scales to the 
effective velocity dispersion (the velocity dispersion within the effective radius) independently of 
the mass or luminosity of the galaxy. The uncertainties in Equation ^ were derived using bootstrap 
resampling of the data points. The cumulative observational uncertainties in this equation are 
0.033 dex in k^, which is substantially smaller than the rms of the fit, implying a substantial 
intrinsic scatter of the "observed mass-to-light ratio" for any given "mass." 

While it is desirable to choose an orthogonal coordinate system which might directly relate 
the observables to underlying physical properties of elliptical galaxies, this conceptualization of 
the FP has a number of problems. First, the quantity R^s is not equivalent to R^, the half-mass 
radius, but instead varies with the observed wavelength. This generally follows from the presence 



of color gradients in elliptical galaxies (e.g., Franx, Illingworth, & Heckman 198£; Peletier et al 



1990a, b), but will be shown explicitly for the case of comparing IZ-band and i('-band effective 
radii in the next paper in this series ( Pahre, de Carvalho, &: Djorgovski 1998| ). It follows that the 



value of Ki, which was intended to create a quantity which is proportional to mass, systematically 
varies as a function of wavelength while mass, of course, does not. Using the terminology of 
Djorgovski et al. (1988), the mapping from R^g to gravitational radius Rg is accomplished in the 
equation R^s = kuRg] the argument here is that the structure function kji is not a constant, as is 
implicitly assumed in the construction of K-space by Bender et al. (1992), but is instead a function 
of wavelength. 

Second, while the evidence is not yet strong, the lowest luminosity ellipticals show no 
detectable color gradients ( peletier et al. 1990aD , suggesting that there could be a dependence of 
the size of the color gradient on luminosity (and hence mass). The size of the color gradients also 
depends on the wavelength sampled (Franx et al. 1989; Peletier et al. 1990a,b). The mappings 
from Ki to mass and K3 to mass-to-light ratio are therefore a function of both wavelength and 
size of the color gradient (which is, in turn, a function of mass). Again using the terminology of 
Djorgovski et al. (1988), this means that the composite structure functions ksR and ksi (see their 
Equations 5-8) are also not constant, as is implied by the construction of K-space, but are instead 
functions of wavelength and mass. 

Third, the use of the central velocity dispersion co in deriving mass at the effective radius 
assumes dynamical homology among elliptical galaxies when mapping do (the central velocity 
dispersion) to aes (the velocity dispersion within the effective radius). Whether or not the 
internal stellar velocity distributions of elliptical galaxies form a homologous family is a point of 
considerable debate. Empirical arguments ( |j0rgensen, Franx, &: Kjaergaard 1995b , Busarello et al. 



the reasons discussed later in this section, we prefer to distinguish between this observed relationship (between 



Ki and K3) and the intrinsic M and (M/L); the intrinsic properties may or may not be fully described by the axes 
Ki and K3. 
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1997) and numerical simulations of dissipation-less merging (Capelato et al. 1995) seem to suggest 
that the way ctq scales to Uefr is a function of galaxy mass or luminosity. The mapping from ki to 
mass and to mass-to-light ratio are therefore a function of mass or luminosity, and possibly a 
function of other physical processes which are currently poorly understood. Once again using the 
terminology of Djorgovski et al. (1988), the mapping from the intrinsic velocity distribution (w^) 
to the observed a in the equation cr^ = kyiv^) is accompanied by a dependence of ky on galaxy 
mass (or luminosity) and systematic deviations from a dynamical homology. 

In summary, because the mapping from the observables to K-space is a function of wavelength, 
luminosity, and deviations from dynamical homology, and furthermore because the mapping 
from K-space to mass and mass-to-light ratio is also a function of wavelength, luminosity, and 
deviations from dynamical homology, we eschew the use of K-space since it is an obfuscation, 
rather than an illumination, of the fundamental physical quantities of elliptical galaxies which 
we wish to understand. At its best, the K-space formalism is merely an intermediate orthogonal 
transformation between the observables [-Reff; (A')eff! ^o] and the desired physical properties 
[M,L,M/L{X)]. 



3.4. The Mg2 o-Q Relation 

The Mg2-cJo relation is a correlation between two distance independent quantities and hence 
useful both as a diagnostic and as a constraint on formation processes for elliptical galaxies as 
a family. Guzman (1995), for example, found that the residuals of the Mg2-(To relation and the 
D~ao relation showed systematic differences between the Hydra-Centaurus region and the Coma 
cluster, thereby suggesting that there are differences between the global properties of galaxies in 
those two environments. 

Only a fraction of galaxies for which do is available also have Mg2 values available. Of 
the entire sample of 301 early-type galaxies in these 16 clusters and groups, only 188 galaxies 
(62%) fit the criteria of logco > 1.8 and have Mg2 measurements. There are six galaxies at 
low (To that show anomalously low Mg2 and are therefore excluded: M32, NGC 3489 in the Leo 
Group, and NGC 4239, NGC 4468, NGC 4476, and NGC 4733 in the Virgo cluster (NGC 4489 
was previously excluded for log do < 1.8). Many of these are dwarf galaxies which are known to 
follow different FP correlations. The criterion used for this exclusion was that all galaxies satisfy 
Mg2 < I (2.2 - logfJo). The best fitting Mg2-ao relation is 

Mg2( mag) = 0.173 log do -0.106 iV = 182 rms = 0.019 mag 

±0.010 ±0.024 ^ ^ 

This relation is plotted in Figure |5[ The slope of this relation is slightly shallower than the value 
of 0.196 ± 0.009 found by J0rgensen (1998), and the scatter is slightly smaller than the 0.025 mag 
of J0rgensen. Inclusion of the six galaxies anomalously low in Mg2 and the one galaxy with 
logcr < 1.8 produces a Mg2-o-o relation Mg2 = 0.188 ± 0.012 log o-q - 0.140 ±0.026 with a scatter of 



0.021 mag; this is closer to, and statistically indistinguishable from, J0rgensen's results. The form 
in Equation |^ will be used, however, as it best represents the properties of the normal elliptical 
galaxies. 



3.5. The Mg2 Near Infrared Fundamental Plane 

An alternative form of the FP was proposed by de Carvalho & Djorgovski (1989): substitute 
a stellar populations indicator, such as the Mg2 index, for the dynamical or mass indicator do in 
the FP relation. The motivation for this is that since Mg2 and <to are strongly correlated with 
each other (as shown above in § p.4D then metallicity could actually be the fundamental physical 
property that causes the slope of the FP to deviate from the virial expectation, de Carvalho & 
Djorgovski also showed that a metallicity sensitive color could be substituted for ao, although that 
approach will not be pursued here due to the heterogeneity of the derived optical-infrared colors 
in Paper II (see Pahre 19981) . p^ The resulting FP in the near-infrared if -band using the Mg2 



index in place of log (Jq is 

logrcfr(arcsec) = 8.3 Mg2+ 0.32A{fiK)eS + Ci iV = 181 rms = 0.172 dex 

±0.9 ±0.015 ^ ' 

The slope of this relation is 8.3 ± 0.9, as predicted by the slope of the Mg2-(To relation combined 
with the standard form of the near-infrared FP, i.e., (1.53 ± 0.08) /(0.173 ± 0.010) = 8.8 ± 0.7. The 
observational uncertainties increase when Mg2 is used instead of log ctq since the latter quantity 
has ~ 3 times greater uncertainty, while the slope has changed by more than a factor of five 
from Equation || to Equation ^. The scatter in the near-infrared Mg2 FP, however, has increased 
by much more than this difference, suggesting that Mg2 is not nearly as good an indicator as 
the velocity dispersion in describing the fundamental and regular physical properties in elliptical 
galaxies that give rise to the FP. The Mg2 index could be identifying real differences in the stellar 
populations among galaxies and hence shows larger scatter when it is substituted into the FP. 

If part of the scatter in the re{f-{fiK)eS~c^o FP can be attributed to differences in stellar 
populations among elliptical galaxies, then the introduction of a stellar populations "correction" 
factor based on the Mg2 index should be able to reduce the scatter of the FP even though a 
small amount of additional observational uncertainty is added in the process. The idea for this 



^ Basically, the fundamental requirement for such an investigation of the "color FP" is to understand the difference 
between an "aperture" color-magnitude relation, the standard form which relies on a color measured in a fixed 
physical aperture size for all galaxies, and a "global" color-magnitude relation, which is evaluated at some fiducial 
scaling radius. Part of the slope of the "aperture" color-magnitude relation is certainly due to the presence of color 
gradients which act in the sense that ellipticals are redder in their centers: the smallest galaxies have their colors 
evaluated at large r/refi where their color is bluer, while the largest galaxies have their colors evaluated at small 
r/vcs where their color is redder. Future work should explicitly distinguish between the two effects of color gradients 
and global color differences in order to place a constraint on the global properties of ellipticals. 
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comes from the attempt by Guzman & Lucey (1993) to construct an "age-independent" distance 
indicator. Here the method will be applied to the near-infrared data. 

Using Bruzual (1983) evolutionary spectral synthesis models, Guzman & Lucey showed that 
the effects of a burst of star formation involving 10% of a galaxy's mass would appear as a change 
of Am/A Mg2 roughly constant for times ^ 1 Gyr or so after the burst. Hence, an offset in 
magnitude Am could be applied to each galaxy independently based on its departure A Mg2 from 
the Mg2-iTo relation. Guzman & Lucey also showed that while the effect in the optical y-band 
was Amv/AMg2 ~ 10, the effect in the near-infrared was much smaller at 

~ 2. (9) 

A Mg2 ^ ^ 

This effect can basically be understood as a filling-in of the Mg2 feature by the addition of a 
continuum flux from hot, young stars. The changes in the Mg2 index are expected to be small, as 
Mg2 is far more sensitive to metallicity than it is to age or IMF ( Mould 1978| ). 



The Guzman & Lucey procedure is repeated here for the Worthey (1994) stellar populations 
models in order to determine if the size of the predicted age effects are similar despite significant 
differences between the Bruzual (1983) and Worthey (1994) models. The model adopted is similar 
to Guzman & Lucey in that it involves 10^ Mq total mass, 90% of which is 15 Gyr old in the 
present day, 10% of which is 5 Gyr old, and all of which has solar metallicity. A similar model 
is investigated that involves 90% of the galaxy being 11 Gyr old in the present day and 10% 
being 5 Gyr old. The results from the Worthey models are plotted in Figure ^ for the UBVRIK 
bandpasses. 

The Worthey models show modest agreement with the Guzman & Lucey calculations based 
on Bruzual (1983) models for the optical bandpasses, although the value of Amag/A Mg2 is 
systematically 30-50% higher. In the /C-band, however, Worthey's models are 3-4 times higher 
in this quantity. The most likely explanation for this is the difficulty that many models have in 
producing enough Mg relative to Fe for very metal-rich populations (see [Worthey, Gonzalez, fc 



Faber 1992| , for example), although more fundamental problems in the treatment of cool stellar 
atmospheres in the infrared for the Worthey models (see [Chariot, Worthey, fc Bressan 1996| ) 
could also be relevant. As a result of this discrepancy, the effects of a late burst of star formation 
involving a small fraction of the galaxy mass on the FP cannot be assumed a priori. 

A better approach is to measure directly the possible contribution of younger stellar 
populations using the FP itself, and then using the observations to constrain the models. 
Combining Equations || and |^ defining tk = , and assuming that the mean surface 

brightness term in the FP is the only one affected by tk yields: 

logreff(arcsec) = a' log ctq + 6' ((Aii^)eff - [ Mg2 - 0.173 log ctq + 0.106]) + q (10) 

where the primed coefficients are the "age-independent" form of the FP. The Guzman &; Lucey 
analysis predicts r/^- ~ 2, while the Worthey models predict tk ~ 7. 
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Unfortunately, a minimization of Equation IC does not reveal an optimal value of rx since 
values of \rK\ > increase the scatter of the equation due to the added measurement uncertainties. 
Put another way, there was no significant improvement in the scatter (by > 2%) of the FP for 
any rx in the range —lOKrx < +10. Since the intrinsic thickness of the near-infrared FP is 
clearly resolved without including a Mg2 term (i.e., Equation ^ and Figure]^, then some of that 
thickness could be due to variations in age among the stellar populations of elliptical galaxies. 
This effect can be viewed from a similar perspective by attempting to correlate the residuals of 
the near-infrared FP with the residuals of the Mg2-o"o relation. This A-A diagram is shown in 
t igure |. This fi gure shows no correlation among the residuals in the direction of the age vectors, 
implying that the resolved intrinsic scatter of both the near-infrared FP and the Mg2-(To relations 
cannot be caused by age effects alone. 

Instead, insight can be gained by looking for changes in the relative distance modulus between 
a given cluster and the Coma cluster. When rx = 2, all 16 clusters and groups show changes in 
their distance moduli relative to Coma of < 0.05 mag, which is smaller than the typical uncertainty 
of 0.1 mag in the distance to a cluster of = 20 galaxies for a scatter of logr^s = 0.085 dex. 
The difference in relative distance modulus between the rx = and rx = 7 cases, however, reach 
as high as 0.12-0.18 mag in the case of several clusters (Virgo, CenSO, and Pegasus), which is 
marginally significant. Since the Worthey (1994) models probably overestimate rx due to their 
difficulty in producing Mg for metal-rich populations, the conclusion is that 0.18 mag is a firm 
upper limit to the effects of age differences on the distance moduli derived using the near-infrared 
FP. This result for the i^-band is similar to that found in the r-band by J0rgensen et al. (1996), 
who found a minimally-significant contribution of r^. = 1.3 it 0.8. They probably found no good 
correlation for the reason given above — that the increase of measurement uncertainties as r 
increases prevents effective minimization during the fitting — hence this may not be a significant 
constraint on a superimposed intermediate age contribution. Repeating the analysis using ry = 10 
for the Abell 2199 and Abell 2634 clusters for the galaxies in the V and K matched catalog 
produces a similarly small change in their distance moduli relative to Coma. In summary, there is 
little evidence that adding a Mg2 term (based on the Mg2-o"o relation) to the near-infrared FP to 
account for age differences in the ellipticals in different clusters causes a significant improvement 
over the zero-point for the relation and hence distances derived from it. 

A similar approach could be pursued by looking at the correlation between the residuals of 
the Mg2-cro relation and the residuals of the k^-ki relation. If differences of Mg2 at a fixed ctq 
indicate differences of stellar populations, and if differences of K3 at a fixed /ti indicate differences 
in {M/L) due to stellar populations effects, then the residuals of these two relations (Equations 10 
and ^) should correlate in a manner that is consistent with stellar populations effects. These 
residuals are plotted against each other in Figure ^. 

This figure is an excellent diagnostic for distinguishing between age and metallicity effects, 
since Mg2 decreases for younger stellar populations while it increases for higher metallicities, 
but (M/Lx) decreases for both. [If the vertical axis were A{M/L) measured at any optical 
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wavelength, then age and metaUicity would instead act nearly parallel.] The lack of any preferred 
correlation along either the age or metallicity vectors in Figure ^ while at the same time having 
a substantial intrinsic scatter for both relations, strongly indicates that both age and metallicity 
variations exist at any given point on both the Mg2 and K3-K1 relations. 

The H/3 index is expected to be a good indicator of the presence of a young stellar component, 
and either the Mg2 or log(Fe) indices should be good indicators of the mean metallicity of 
the stellar content. The large intrinsic scatter between the H/3 and Mg2 indices, the lack of 
a correlation altogether between H/3 and log(Fe), the strong correlation between Mg2 and uo, 
and the weak correlation between log(Fe) and (Tq ( lJ0rgensen 199§| ) all indicate that there exist 
significant variations in both age and metallicity for any given value of fio. This is fully consistent 
with the argument above based on the residuals of the Mg2-f7o and relations. 



3.6. The Faber Jackson Relation 

The correlation between luminosity and central velocity dispersion for elliptical galaxies was 
first noticed by Faber &: Jackson (1976). If we fit the relation L oc CTq, this is equivalent to fitting 
Ktot = —2.5a log (To + b. The best fitting Faber-Jackson relation for the iiT-band data is 

Mk = Ktot - 34.91 + 5 log /i75 = -10.35 log uq N = 252 rms = 0.93 mag 

±0.55 ^ '* 

assuming Hq = 75 km s~^ and that the Coma cluster {cz = 7200 km s~^) is at rest with respect to 
the Hubble flow. The relation is plotted in Figure The scatter of this relation is significantly 
smaller in the Coma cluster alone (rms= 0.72 mag). 



3.7. The Modified Faber-Jackson Form of the FP 

Since there is substantial scatter in the Faber-Jackson relation due to variations in surface 
brightness among galaxies at a given luminosity and central velocity dispersion, an alternate form 
of the FP is to substitute K^ot for reff. This will be referred to as the "modified Faber-Jackson" 
relation, as it adds the additional (//)efr term to the Faber-Jackson relation. The form of this 
equation is 

Ktot = a'log(To(km s"^) + 6'(^)efr(mag arcsec"^) + c- (12) 

where the primed coefficients are used here for the modified Faber-Jackson relation. Since in 
the case of a pure de Vaucouleurs profile irtot(mag) = — 51ogreff + (^A')efr — 1.995, the value 
logTeff = a log (To + b{fiK)eS + Constant can be substituted for r^g resulting in a = —5a' and 
b = 0.2(1 — b'), thereby relating Equations || and 12 to each other. The best-fitting relation of the 
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modified Faber-Jackson form of the FP is: 

Mx + 51og/i75 = -8.16 log do- 0.585(/Xi^)eff + Ci iV = 251 rms = 0.51 mag 

±0.47 ±0.062 ^ ' 

which has a scatter only 10% larger than the standard form of the near-infrared FP given in 
Equation |2|. This equation represents the scaling relation Lk oc a^'"^^^^'^^ {T,k)^q^'^^^'^^ ■, which is 
fully equivalent within the uncertainties to the standard form of the near-infrared FP. Since the 
uncertainty on Ktot ± 0.60(//x)eff is 0.068 mag, the total observational uncertainties in Equation § 



and M are similar. Hence Equation 13 also shows substantial intrinsic scatter in the properties of 



elliptical galaxies at any point along the FP. 



3.8. The Kormendy Relation 

Effective radius and mean surface brightness, which are two of the three terms in the FP, are 
correlated with each other ( Kormendy 1977| ). The best fitting Kormendy relation for the X-band 



data is 

logi?eff/ire = 0.244(/Xi^)eff - 3.637 = 269 rms = 0.227 dex 

±0.029 ^ ' 

assuming Hq = 75 km s^^ and that the Coma cluster [cz = 7200 km s~^) is at rest with respect to 
the Hubble flow. The relation is plotted in Figure The scatter of this relation is significantly 
smaller in the Coma cluster alone (rms= 0.198 dex). Measurement errors in r^^ and (/i)eff are 
correlated and act in a direction nearly parallel to the Kormendy relation, but are not nearly 
large enough to account for the spread in galaxy properties along the relation. Changes in 



luminosity are skewed with respect to this relation and shown in Figure |lj. For this reason it is 
necessary that a magnitude-limited sample be defined in a consistent way for all clusters studied 
before conclusions based on changes in the zero-point (due to distances for nearby clusters or 
evolutionary brightening for higher redshifts) can be made. 



3.9. The Radius Luminosity Relation 

The correlation between the effective radius and total magnitude for elliptical galaxies has 



long been used for distance scale work and especially cosmological tests (see Sandage k, Perelmuter 



1990 and references therein). The best fitting radius-luminosity relation for the ii'-band data is 
± 5 log /i75 = -4.40i2eff - 22.31 N = 269 rms = 0.88 mag 

(15) 

±0.26 ^ ' 

assuming Hq = 75 km s^^ and that the Coma cluster {cz = 7200 km s~^) is at rest with respect 
to the Hubble flow. The relation is plotted as Figure |l^. There is intrinsic scatter to this relation 
that is a result of the variation in surface brightness at a given radius and luminosity; lines of 
constant surface brightness are plotted in the figure to demonstrate this effect. 
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4. Exploring Simple Models for the Origins of the Elliptical Galaxy Scaling 

Relations in the Near Infrared 



The near-infrared FP has been shown in ^3.1| to be represented by the scahng relation 
Teff oc a^^^^^'^^ {T,k)cS "^^^^'^^ ■ This relation shows a significant deviation from the optical 
forms of the FP: r^s oc o-^-24±0-07^5]^-^o.82±o.02 (|j0rgensen, Franx, fc Kjaergaard 199^) and 
reff oc ct1-38±o.04^5.^-^o.82±o.03 ( [Hudson 6^ d. 1997D in the i?-band; or r^s oc ay^J:)J-'^^ in the 
l^-band m uzman, Lucey, Sz Bower 19931) . There are two simple conclusions to draw from these 
data: (1) the slope of the near-infrared FP deviates from the virial expectation of r^s oc c7o(S)^^, 
and (2) the slope of the FP increases with wavelength. A third insight derives from the fact that 
the scatter of the FP is very similar at all wavelengths. These three points are sufficient to discuss 
several simple models for the physical origins of the FP. 

The age-met allicity model of Worthey, Trager, & Faber (1995) — based on the form of the FP 
in the optical, various line indices, and simple stellar populations model comparisons — incorrectly 
predicts that the near-infrared FP should follow the virial form. Another model, that the FP 
slope is caused by deviations of the velocity distributions of elliptical galaxies from a homologous 
scaling family (Capelato et al. 1995), cannot account for the variations of the slope of the FP 
with wavelength. If this breaking of homology has its origin in dissipation-less merging, then this 
effect also cannot explain the correlation between Mg2 and ctq. A final model, which suggests that 
deviations of the light distributions of elliptical galaxies from the de Vaucouleurs r^/^ form is the 



cause of the FP slope, is unable to account for the slope of the FP in the optical (Graham & 



Colless 1997), and for the same reasons it cannot explain the near-infrared FP slope. 



The deviation of the slope of the near-infrared FP from the virial expectation, assuming 
homology and constant M/L among ellipticals, is a very significant result. This requires a 
breakdown of one or both assumptions: either M/L is systematically varying along the FP, or 
elliptical galaxies are systematically deviating from a homologous scaling family. If age is the 
stellar populations parameter which causes variations in the slope of the FP with wavelength, 
then age alone might possibly produce the slope of the K-hand FP. This conclusion, however, is 
severely limited by the possibility of homology breaking along the elliptical galaxy sequence. 

Allowing for structural deviations from homology, in the form of a Sersic r^/" profile, does not 
appear to cause significant changes to the slope of the FP for high S/N, y-band data in the Virgo 
cluster ( praham fc Colless 1997 ). Instead, allowing for dynamical deviations from homology, via 



galaxy to galaxy variations in the mapping from do to cJefr, appears to cause significant changes 
in the slope of the FP ( [Busarello et al. 1997 ; cf. Graham & Colless 1997). Busarello et al. 



found a relationship between the velocity dispersions to be logdo = (1.28 ± 0.11) log cJeff — 0.58. 
Substituting for log ctq into the i('-band FP solution in Equation ^ produces 

log reff = 1.96 log creff+ 0.314(/ix)eff + Constant 

±0.20 ±0.011 ^ ' 

which is statistically consistent with the virial expectation of reg oc cr^. This argument suggests 
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that the deviation of the near-infrared FP from the virial expectation can be fully explained by 
systematic deviations of the velocity structure of elliptical galaxies from a homologous family, 
removing the requirement of large age spreads among elliptical galaxies. 

Since either dynamical non-homology or large age spreads could produce the slope of the 
near-infrared FP, it is impossible to distinguish between these two simple models without further 
analysis. In addition, any model which incorporates either age or dynamical deviations from a 
homology along with metallicity variations cannot be excluded in this simple analysis, either. This 
strongly suggests that a much more detailed analysis, along with a more complicated model with 
several different variables, is necessary to explain the global properties of elliptical galaxies. 

One last ad hoc model can be constructed in which there is a conspiracy between metallicity 
and age effects that act in a manner to keep the FP thin. In this model, there can be a large spread 
in age and metallicity at any given point of the FP — under the constraint that the two effects of 
age and metallicity work opposite to each other and thereby cancel out to maintain a small scatter. 
While this model would work at optical wavelengths, the independence of near-infrared light to 
metallicity would cause the thinness of the optical FP to break down into a thick, near-infrared 
FP. Since the near-infrared FP has similar observed and intrinsic thicknesses when compared 
to the optical FP, especially when considering the additional observational uncertainties on r^s 
caused by ctq for the steeper slope of the near-infrared FP, this model can be excluded. 



5. Discussion 

The near-infrared FP that has been constructed in this paper has several important 
properties: (1) it deviates from the virial expectation (assuming constant M/L and homology); (2) 
it is steeper than the optical FP relations; (3) it has a similarly small scatter when compared to 
the optical FP relations; and (4) it has a small, but significantly resolved, intrinsic scatter. These 
observational constraints are sufficient to exclude a number of simple models for the origin of the 
FP, but they do not provide unique discrimination between composite models which include either 
age, systematic deviations from dynamical homology, or both. Small additional contributions due 
to metallicity variations are also possible. 

Better insight is gained by including the Mg2 index into the analysis. The Mg2 form of the 
FP has much larger scatter than the standard form, which argues that Mg2 does not uniquely 
specify the depth of the potential well for each galaxy. This is entirely consistent with the resolved 
intrinsic scatter of the Mg2-o"o relation. If some physical process like galactic winds (Yoshii & 



Arimoto 1987] ) caused the metallicity and potential well for all elliptical galaxies to behave like 



a one-parameter family, then some other physical property, such as dissipation-less mergers or a 
large scatter in formation times, would be required to produce the small intrinsic scatter of the 
Mg2-o"o relation and the large intrinsic scatter of the Mg2 form of the FP. 



The near-infrared FP has the unique property that the X-band light is virtually independent 
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of metallicity. For this reason, residuals of the if -band FP (or the if -band relationship between 

Ki and K3) and the Mg2-(To relations could provide a strong discrimination between age and 
metallicity effects. The lack of any clear correlation between these residuals implies that neither 
age nor metallicity is a unique contributor to the intrinsic scatter of the FP or the Mg2-(To 
relations. 

The Fundamental Plane is not just a simple correlation of the observed properties of elliptical 
galaxies, but rather a unique tool for studying the intrinsic physical properties spanned by these 
galaxies. The remarkable homogeneity of properties of elliptical galaxies that is implied by the 
regularity and thinness of the optical FP is clearly reproduced by their similarly regular properties 
in the near-infrared bandpass. 
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Table 1. Fits for Each Cluster or Group for the Near-Infrared FP 



Cluster or Simultaneous Fit Individual Fits Constrained b = 0.314 Fits 

Group Ci N rms a Aa b Ab rms a Aa rms 

(dex) (dex) (dex) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Coma 


-7, 


,950 


60 


0, 


,086 


1, 


.33 


0. 


19 


0, 


.302 


0. 


,03 


0. 


,082 


1, 


.57 


0, 


.15 


0, 


.088 


A 194 


-7, 


,734 


16 


0, 


107 


1, 


.57 


0. 


,21 


0, 


.254 


0. 


,05 


0. 


106 


1, 


.60 


0, 


.16 


0, 


.110 


A2199 


-8, 


,128 


23 


0. 


,093 


1, 


.53 


0. 


,22 


0, 


.342 


0. 


,03 


0. 


,088 


1 


.40 


0, 


.16 


0, 


.086 


A2634 


-8, 


,028 


15 


0, 


,076 


1, 


.19 


0. 


,74 


0, 


.292 


0. 


,06 


0. 


,061 


1, 


.24 


0, 


.30 


0, 


.063 


Cen45 


-7, 


,543 


6 


0, 


,071 


































CenSO 


-7, 


,526 


14 


0, 


124 


1, 


.72 


0. 


,41 


0, 


.299 


0. 


,08 


0. 


123 


2 


.05 


0, 


.39 


0, 


.139 


Fornax 


-7, 


,274 


15 


0, 


137 


2, 


.56 


0. 


,65 


0, 


.339 


0. 


,06 


0. 


156 


2 


.11 


0, 


.30 


0, 


.128 


Hydra 


-7, 


,669 


17 


0, 


,086 


1, 


.76 


0. 


,34 


0, 


.344 


0. 


,03 


0. 


,080 


1 


.77 


0, 


.15 


0, 


.086 


Klemola 44 


-8, 


,041 


11 


0, 


,067 


1, 


.50 


0. 


,58 


0, 


.309 


0. 


,05 


0. 


,068 


1 


.74 


0, 


.28 


0, 


.069 


Pegasus 


-7. 


,580 


4 


0. 


,048 


































Perseus 


-7, 


,802 


19 


0, 


100 


1, 


.98 


0. 


,59 


0, 


.310 


0. 


,05 


0. 


125 


1, 


.66 


0, 


.33 





.104 


Pisces 


-7, 


,723 


11 


0, 


,087 


1, 


.04 


0. 


,27 


0, 


.350 


0. 


,05 


0. 


,055 


1 


.17 


0, 


.19 


0, 


.055 


Virgo 


-7, 


175 


27 


0. 


,115 


1, 


.77 


0. 


,25 


0, 


.374 


0. 


,03 


0. 


120 


1 


.62 


0, 


.12 


0, 


.118 


Eridanus 


-7, 


,312 


5 


0, 


,061 


































Leo 


-6, 


,932 


3 


0, 


,253 


































N5846grp 


-7, 


,436 


5 


0. 


101 



































Note. — (1) The FP fits in this table are to the form of Equation gj. (2) The simultaneous fit for columns 
2-4 corresponds to the solution in Equation ^, allowing only the intercepts Ci to vary between clusters. (3) The 
individual cluster FP fits in columns 5-9 are for only those 11 clusters with numbers of galaxies TV > 10. (4) The 
constrained, individual cluster fits of columns 10-12 were obtained by fixing b = 0.314. (5) The rms in all cases is 
evaluated along the log r^ff axis. 



-20- 



REFERENCES 

Baum, W. A. 1959, PASP, 71, 106 

Bender, R., Burstein, D., & Faber, S. M. 1992, ApJ, 399, 462 
Bower, R. G., Lucey, J. R., & Ellis, R. S. 1992, MNRAS, 254, 601 
Bruzual, A. G. 1983, ApJ, 273, 105 

Busarello, G., Capaccioli, M., Capozziello, S., Longo, G., & Puddu, E. 1997, A&A, 320, 415 
Capelato, H. V., de Carvalho, R. R., & Carlberg, R. G. 1995, ApJ, 451, 525 

Capelato, H. V., de Carvalho, R. R., & Carlberg, R. G. 1997, in Galaxy Scaling Relations: Origins, 
Evolution, and Applications, Proceedings of the Third ESO-VLT Workshop, eds. L. N. da 
Costa & A. Renzini (Springer-Verlag: Berlin), 331 

de Carvalho, R. R., & Djorgovski, S. 1989, ApJ, 341, L37 

Chariot, S., Worthey, G., & Bressan, A. 1996, ApJ, 457, 625 

Ciotti, L., Lanzoni, B., & Renzini, A. 1996, MNRAS, 282, 1 

Djorgovski, S. G., de Carvalho, R., & Han, M.-S. 1988, in The Extragalactic Distance Scale, ASP 
Conf. Ser. Vol. 24, eds. S. van den Bergh & C. J. Pritchet (San Francisco, ASP), 329 

Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 

Djorgovski, S., &: Santiago, B. X. 1993, in Proceedings of the ESO/EIPC Workshop on Structure, 
Dynamics, and Chemical Evolution of Early-Type Galaxies, ed. J. Danziger, et al., ESO 
publication No. 45, 59 

Dressier, A. 1984, ApJ, 281, 512 

Dressier, A., Lynden-Bell, D., Burstein, D., Davies, R. L., Faber, S. M., Terlevich, R. J., &; 
Wegner, G. 1987, ApJ, 313, 42 

Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 

Franx, M., Illingworth, G., & Heckman, T. 1989, AJ, 98, 538 

Graham, A., & Colless, M. 1997, MNRAS, 287, 221 

Guzman, R., Lucey, J. R., & Bower, R. G. 1993, MNRAS, 265, 731 

Guzman, R. 1995, in Proceedings of the Heron Island Workshop on Peculiar Velocities in the 
Universe, http : //qso . lanl . gov/ "heron/ 



- 21 - 



Hudson, M. J., Lucey, J. R., Smith, R. J., & Steel, J. 1997, MNRAS, in press 

JefFerys, W. H., Fitzpatrick, M. J., McArthur, B. E., and McCartney, J. E. 1987, "GaussFit: A 
System for Least Squares and Robust Estimation," The University of Texas at Austin 

J0rgensen, I. 1998, MNRAS, in press 

J0rgensen, L, Pranx, M., k Kjaergaard, R 1995a, MNRAS, 273, 1097 
J0rgensen, L, Pranx, M., & Kjaergaard, P. 1995b, MNRAS, 276, 1341 
J0rgensen, I., Franx, M., & Kjaergaard, P. 1996, MNRAS, 280, 167 
Kormendy, J. 1977, ApJ, 218, 333 
Lucey, J. R., & Carter, D. 1988, MNRAS, 235, 1177 

Lynden-Bell, D., Faber, S. M., Burstein, D., Davies, R. L., Dressier, A., Terlevich, R. J., &: 
Wegner, G. 1988, ApJ, 326, 19 

Mould, J. R. 1978, ApJ, 220, 434 

Murphy, D. C, Pcrsson, S. E., Pahre, M. A., Sivaramakrishnan, A., & Djorgovski, S. G. 1995, 
PASP, 107, 1234 

Pahre, M. A. 1998, ApJS, submitted 

Pahre, M. A., de Carvalho, R. R., & Djorgovski, S. G. 1998, AJ, submitted 

Pahre, M. A., & Djorgovski, S. G. 1997, in The Nature of Elliptical Galaxies, Proceedings of the 
Second Stromlo Symposium, eds. M. Arnaboldi, G. S. Da Costa, &l P. Saha, ASP Conf. 
Ser. Vol. 116, (San Francisco: ASP), 154 

Pahre, M. A., Djorgovski, S. G., & de Carvalho, R. R. 1995, ApJ, 453, L17 

Pahre, M. A., Djorgovski, S. G., &; de Carvalho, R. R. 1997, in Galaxy Scaling Relations: Origins, 
Evolution, and Applications, Proceedings of the Third ESO-VLT Workshop, eds. L. N. da 
Costa k. A. Renzini (Springer- Verlag: Berlin), 197 

Peletier, R. F., Davies, R. L., Illingworth, G. D., Davis, L. E., & Cawson, M. 1990a, AJ, 100, 1091 

Peletier, R. F., Valentijn, E. A., & Jameson, R. F. 1990b, A&A, 233, 62 

Persson, S. E., Frogel, J. A., & Aaronson, M. 1979, ApJS, 39, 61 

Recillas-Cruz, E., Carrasco, L., Serrano, P. G., & Cruz- Gonzalez, I. 1990, A&A, 229, 64 
Recillas-Cruz, E., Carrasco, L., Serrano, P. G., & Cruz-Gonzalez, L 1991, A&A, 249, 312 



- 22 - 



Sandage, A. 1972, ApJ, 176, 21 

Sandage, A., & Visvanathan, N. 1978, ApJ, 223, 707 
Sandage, A., & Visvanathan, N. 1978, ApJ, 228, 81 
Sandage, A., & Perelmuter, J.-M. 1990, ApJ, 361, 1 

Smith, R. J., Lucey, J. R., Hudson, M. J., & Steel, J. 1997, MNRAS, in press 
Terlevich, R. J., Davies, R. L., Faber, S. M., & Burstein, D. 1981, MNRAS, 196, 381 
Tonry, J. L., & Davis, J. 1981, ApJ, 246, 680 

Tonry, J. L., Blakeslee, J. P., Ajhar, E. A., & Dressier, A. 1997, ApJ, 475, 399 
Visvanathan, N., & Sandage, A. 1977, ApJ, 216, 214 
Worthey, G. 1994, ApJS, 95, 107 

Worthey, G., Gonzalez, J. J., & Faber, S. M. 1992, ApJS, 398, 69 

Worthey, G., Trager, S. C, & Faber, S. M. 1995, in Fresh Views of EUiptical Galaxies, ASP Conf. 
Ser. Vol. 86, eds. A. Buzzoni &: A. Rcnzini (San Francisco: ASP) 203 

Yoshii, Y., & Arimoto, N. 1987, A&A, 188, 13 



This preprint was prepared with the AAS lAT^^X macros v4.0. 



- 23 - 




456456456456 



log ao(km/s) + 0.206 </^K>eff 

Fig. 1. — The Fundamental Plane in the near-infrared for the 16 clusters and groups in the 
simultaneous fit represented by the solution of Equation |2| and the intercepts in column 2 of Table |l|. 
The FP is described by the scaling relation r^.^ oc aQ-^^ {Y,k)cs "^^ with a scatter of 0.096 dex in 
logTeff; the scatter is reduced by 10% of the galaxies with (Tq < 100 km s~^ are excluded. The 
fitted galaxies are plotted as solid symbols, while those excluded from the fit (log ao < 1.8, late-type 
morphology, r^fi < 1 arcsec, or background/foreground in the Virgo cluster) are plotted as open 
symbols. The FP fit is plotted in each panel as a solid line. 
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Fig. 2. — (a) The Fundamental Plane in the ET-band for the combined 16 cluster and group sample, 
seen edge-on along its long side. The symbols are as in Figure |l[ The ordinate is in units of kpc 
assuming Hq = 75 km s Mpc~^. (b) The FP in the i^-band for the combined 16 cluster and group 
sample, seen edge-on along its short side. In this view of the FP, the observationally-correlated 
measurement errors in i?cff and {^K)eS (ordinate) are separated from the independent measurement 
uncertainties in gq (abscissa); the typical measurement uncertainties are shown in the lower right- 
hand corner of the panel, (c) The FP seen face-on. Galaxies do not uniformly populate this planar 
surface. While the i^-band data in this paper are not drawn from a strictly magnitude-limited 
sample, they do behave as though a K^ot ^ 13 mag limit (dotted line) were imposed. Most of the 
galaxies are found to have 15 < {fiK)eS < 18 mag arcsec"^ (long-dashed lines), although there are 
no clear selection effects causing this distribution of galaxy properties. Furthermore, there are no 
galaxies with properties in the upper-right portion of the figure, which could be caused by the lack 
of galaxies with velocity dispersions fio > 400 km s"^ (short-dashed line). The vectors drawn in 
the lower left-hand corner show the direction in which each of the observed quantities varies along 
the plane. 
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Fig. 3. — The Dk-ctq relation for the galaxies in the 16 clusters and groups of the survey, plotted 
in physical units. The scatter for this relation is 0.112 dex in log Dk, which is slightly higher than 
the FP itself, since the Dk-ctq is not quite an edge-on view of the FP. 
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Fig. 4. — The i^-band FP viewed in the K-space perspective. This coordinate system (defined 
by Bender, Burstein, fc Faber 1992| ) is given in Equation |5|, and was designed such that ki is 
roughly proportional to the logarithm of mass and K3 is roughly proportional to the logarithm 
of mass-to-light ratio. The fit in the top panel corresponds to the "observed" scaling relations 
(M/Lk) oc M°-15=^°-°^ or {M/Lk) oc under the assumptions that there are no color 

gradients in elliptical galaxies and dynamical homology is preserved within the family of elliptical 
galaxies. 
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Fig. 5. — The Mg2-cro relation for the 182 galaxies with Mg2 measurements. The scatter about this 
relation is 0.019 mag in Mg2, which is significantly larger than the typical measurement uncertainties 
of 0.013 mag (shown in lower right-hand corner of figure). 
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Fig. 6. — The near-infrared FP with the Mg2 index substituted for the velocity dispersion. 
This figure is plotted to the same scale as Figure |2|(a), hence a direct comparison of these two 
figures demonstrates how the scatter of the FP relation has increased by a factor of two by the 
substitution of Mg2 for ctq. Only a small part of this increase in scatter can be attributed to the 
larger measurement uncertainties of Mg2 compared to (Tq, hence the correlation plotted in this 
figure cannot be an edge-on view of the FP. 
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Fig. 7. — The ratio of the change in magnitude Am to the change in the Mg2 index for the 

UBVRIK bandpasses using the models of Worthey (1994). The two cases considered both have a 
90% (by mass) old stellar population component of 11 Gyr (a) and 15 Gyr (b) and [Fe/H]= in 
addition to a 10% young stellar population that is 5 Gyr old at the present day. 
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Fig. 8. — The residuals of the near-infrared FP plotted against the residuals of the Mg2-(To relation. 
For the purpose of this comparison, only those galaxies with good Mg2 measurements (see § |3.4[) 
were used to fit the near-infrared FP; the FP fit for these 182 galaxies is rgg oc a^'^^ {'Sk)^^''^^ ■ 
A correlation due to the measurement errors of log (Tq would act in the direction of the vectors in 
the lower-right of the panel labeled Actq. If a galaxy had a starburst involving 10% of its mass at 
5 Gyr before the present day, this would produce an offset as shown for the Bruzual (1983; B83) 
and Worthey (1994; W94) models. Since there is no such correlation along these model vectors 
between the two sets of residuals, this implies that there is no age effect, as traced by the Mg2-o"o 
relation, which is the sole cause of the intrinsic scatter of the near-infrared FP. 
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Fig. 9. — The residuals from the Mg2-cTo relation (Equation [lO|) plotted against the residuals from 
the K3-K1 relation (Equation ^). Both relations show significant intrinsic scatter which could be due 
to variations in stellar populations at any given point on either relation. As shown by the vectors 
in the upper part of the figure, the Bruzual & Chariot (1996) models show that this diagram is a 
powerful diagnostic for separating the effects of age and metallicity. These models show the effect 
for changing the age from 11 to 8 Gyr for a solar metallicity population, and for changing [Fe/H] 
from to 0.25 dex in a 11 Gyr old population. The effect of correlated errors in ao are shown in the 
lower-right of the figure. Variations in both age and metallicity at any given mass or luminosity 
appear to be necessary to explain the intrinsic scatter in this diagram and the lack of correlation 
along either vector in the upper right of the figure. 
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Fig. 10. — The Faber-Jackson relation between luminosity and central velocity dispersion. The 
best-fitting relation is Lk oc ctq'^^^'''^^ with a large scatter of 0.93 mag. The scatter is significantly 
smaller in the Coma cluster at 0.72 mag. 
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Fig. 11. — The "modified Faber-Jackson" form of the FP. The best-fitting relation is Lk oc 
^3.26±o.i9^5^^^-o.59±o.06 ^-^^ ^ scatter of 0.51 mag. This form of the FP is nearly identical to 
that of Figure |2| and Equation although it shows 10% larger scatter primarily due to larger 
observational uncertainties. 
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Fig. 12. — The Kormcndy relation between effective radius and mean surface brightness. The best- 
fitting relation is Res oc with a scatter of 0.227 dex in logi?efr, which is 2.4 times 
worse than the scatter of the FP (which has the additional cto term). The scatter is significantly 
smaller in the Coma cluster at 0.198 dex. Measurement errors in R^s and {n)eS are correlated nearly 
along the relation, as shown by the representative error bars in the lower-left of the figure. The 
measurement uncertainties perpendicular to the relation are only 0.015 dex in log i?cfr ~ 0.32 (;U;^)cfr . 
Luminosity changes act perpendicular to the relation, and representative lines of Mk = —22 mag 
to Mk = —27 mag are shown as dotted lines. 
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Fig. 13. — The relation between i^-band total luminosity and effective radius given by Lk esc 

-^eff^^^^*^- There is substantial intrinsic scatter in this relation due to variations in surface 

brightness. The dotted lines each represent the variation of Mk and R^s at constant surface 
brightness. 



